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ABSTRACT 

Context. There is a lively debate about the solar oxygen abundance and the role of 3D models in its recent downward revision. These 
models have been tested using high resolution solar atlases of flux and disk-centre intensity. Further testing can be done using centre- 
to-limb variations. 

Aims. Using high resolution and high S/N observations of neutral oxygen lines across the solar surface we seek to test if the 3D 
and ID models reproduce their observed centre-to-limb variation. In particular we seek to assess whether the latest generation of 3D 
hydrodynamical solar model atmospheres and NLTE line formation calculations are appropriate to derive the solar oxygen abundance. 
Mettiods. We use our recent observations of Oi 777 nm, Oi 615.81 nm, [Oi] 630.03 nm and nine lines of other elements for five 
viewing angles 0.2 < ji < I of the quiet solar disk. We compare them with the predicted line profiles from the 3D and ID models 
computed using the most up-to-date line formation codes, line data and allowing for departures of UTE. The centre-to-limb variation 
of the O I 777 nm lines is also used to obtain an empirical correction for the poorly known efficiency of the inelastic collisions with 
Hi. 

Results. The 3D model generally reproduces the centre-to-limb observations of the lines very well, particularly the oxygen lines. 
From the O i 777 nm lines we find that the classical Drawin recipe slightly overestimates H i collisions (we find the best agreement 
with 5 H ~ 0.85). The limb observations of the Oi 615.82 nm line allow us to identify a previously unknown contribution of molecules 
for this line, prevalent at the solar limb. A detailed treatment of the [O i] 630.03 nm line including the recent nickel abundance shows 
that the 3D modeling provides an excellent agreement with the observations. The derived oxygen abundances with the 3D model are 
8.68 (777 nm lines), 8.66 (630.03 nm line) and 8.62 (615.82 nm fine). 

Conclusions. These additional tests have reinforced the trustworthiness of the 3D model and line formation for abundance analyses. 
Key words, line: formation - Sun: photosphere - Sun: granulation ~ Sun: abundances - Convection 



1. Introduction 

Oxygen is arguably one of the most important elements in the 
Universe. After hydrogen and helium it is the most common el- 
ement, projecting its importance to many fields in astrophysics: 
from stellar abundances to inter stellar medium and galactic evo- 
lution. The solar oxygen abundance is often used as a reference 
in many studies. Being a volatile element, its meteoritic abun- 
dance is not representative. Hence the most reliable probe for 
the reference oxygen abundance is the solar photosphere. 

With high-quality data readily available, a long standing ex- 
pertise in solar observations and much improved knowledge of 
stellar atmospheres in the last decades, one would expect that 
a measurement as fundamental as the solar photospheric abun- 
dance of oxygen would be firmly established. But, on the con- 
trary, the solar oxygen abundance has been hotly debated in re- 
cent times. From the high (and widely used) value of log sq - 
8.93 + 0.04 ( [Anders & Grevesse|[T9g9l l, 
abundance has been revised downward 
log eo - 8.66+0.05 ( Asplund et al. 2004 ). The pivotal causes for 
this change have been the proper treatment of the statistical inho- 
mogeneities caused by the solar granulation (making use of a 3D 



the proposed oxygen 
to the low value of 



photosphere model); the proper treatment of departures from lo- 
cal thermodynamical equilibrium (LTE) and improved line data 
(including better blend identification). A low oxygen abundance 
has caused significant grievances among the solar/stellar interior 
modelling community by undoing the almost perfect agreement 
between the solar interior models and helioseismology, which 
has sparked much debate. Over the last years the downward re- 
vision of the solar oxygen has been supported by some studies 



(e.g.jSocas-Navarro & Norton 2007 


Melendez & Asplund 2008| 


but also contested by others (e.g. 


Centeno & Socas-Navarro] 



2008 Ayres|[2008 1. Some recent studies find intermediate oxy 



gen abundances (e.g. Caffau et al.|2008[ l. 



A measure of the photospheric oxygen abundance is made 
difficult for several reasons. Of the few atomic oxygen lines 
available in the solar spectrum, some are very weak, others are 
significantly blended. An example is the popular [O i] 630.03 nm 
line, which is weak and has a non-negligible blend with a nickel 
line. On the other hand, departures from LTE are also important 
for some O i lines - the O i 777 nm triplet lines being the typical 
example. These lines are known to show significant departures 
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from LTE, requiring the use of detailed NLTE line formation that 
requires detailed atomic input data which may not exist. 

For the NLTE modeling the collisions with neutral hydro- 
gen are particularly relevant. There is disagreement in the liter- 
ature regarding the importance of inelastic collisions with Hi. 
At this point, the lack of experimental data or quantum me- 
chanical calculations makes its estimation complicated. A cus- 
tomary approach is to adopt a generalization of the classical 
Drawin formula ( Drawin|[l968[ l, using the recipe of [Steenbock 



O I 615.81 nm 



|& Holweger| (|1984[ ), often scaled by an empirical factor Sh\ 
note that |Lambertr ( |1993| l have corrected a mistake in these 
original formulae. However, there are different views regarding 
which scaling factor to use for oxygen. |Nissen et al._(2002j and 



Asplund et al. ( 2004[ l chose Sh = 0, based on evidence that 



the Drawin formula overestimates the Hi collision efficiency 
with some atoms for which experimental or detailed quantum 
mechanical data exists (e.g. |Barklem et al.||2003] l. |Caffau et al. 
( |200 8 ), on the other hand, have adopted a seemingly ad-hoc 
value of Sh - 1/3. 

For the solar Oi 777 nm lines, the adopted recipe for Hi 
collisions has a significant effect on the line shape and strength 
and consequently, on the derived oxygen abundance. Using the 
same 3D model of Asplund et al. ( 2004[ l, AUende Prieto et al. 
( 2004| l empirically found that 5h = 1 agreed somewhat better 
than 5h = while the LTE case could be ruled out. 

In addition to the line formation physics, the centre-to-limb 
variation of the lines provides a robust test of the model atmo- 
spheres by probing the depth variation of the source functions, 
as seen from the Eddington-Barbier approximation. In this work 
we study atomic oxygen lines using observations of their centre- 
to-limb variation. We study the effects of using different atmo- 
spheric models, departures from LTE and properly accounting 
for blends. Our aim is to provide additional observational tests 
of some of the models used to infer the oxygen abundance, using 
new solar observations. 

We have obtained high spatial and spectral resolution ob- 
servations of oxygen and other lines across the solar surface. 
Similar observations of oxygen hnes have been obtained in the 
past (notabl y Muller et al.||1968[ |AItrock|[T968{ [Allende Pri"eto 



|et al.||2004 i. However, available data from early works is lim- 



ited to equivalent widths and the more recent work of Allende 



Prieto et al. 



( 2004 1 covers only the O i 777 nm triplet hnes. 

Furthermore, the long slit used by Allende Prieto et aT] ( |2004| l 
increases the uncertainty of the fi value for their limb data. The 
observations detailed in the present work include five neutral 
oxygen lines: the Oi 615.81 nm, [Oi] 630.03 nm and the three 
1 777 nm over five positions in the solar disk. In addition, lines 
from other elements are also included in the observed spectral 
regions. 

This paper is organized in the following manner. Our ob- 
servations are briefly outlined in the next section. In Section [3] 
we detail the model atmospheres and the line formation codes 
used, which are compared with the observations in Section |4] 
Conclusions are made in Section |5] 



2. Observations 

2.1. Overview 



We make use of the observations of |Pereira et al.| ( |2009l hereafter 
[Paper I[ ), which we refer to for a detailed description of the pro- 
gramme, instruments and reduction. Whereas Paper I is focused 
on high spatial resolution, here we study the centre-to-limb vari- 
ations of the spatial and temporaly mean spectra. Thus, in terms 
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Fig. 1. Observed profiles at disk-centre and limb. For the weak 
615.81 and 630.03 nm lines the location is indicated. In this fig- 
ure the wavelengths for the limb spectra have been adjusted to 
compensate for solar rotation and other systematics. 



of the observations the main difference is that we average the 
spectrograms in space and time to obtain a mean spectrum for 
each position in the solar disk. 

Using 25 spectrograms for each /i 1 position and 50 spec- 
trograms for i-i - I means that our mean spectra comprise an 
average of more than 20000 spectra (> 40000 for disk-centre). 
The S/N ratio for the spectrum of one spatially averaged spec- 
trogram is about 700. The total S/N would be V25 times that 
value, if the images were all independent - which is not the case 
because the time separation does not always guarantee they sam- 
ple a different granulation pattern. In any case, a conservative 
estimate of the S/N should be at least 1500 - more than enough 



for the present analysis. Unlike in Paper I no Fourier filtering 



is applied, since most of the photon noise is eliminated by the 
averaging of many spectra. 

It is worth to note that for the jj. x 0.4 set some facular fea- 
tures were visible on the slit-jaw images, though the slit did not 
cross them. Because we want to study quiet regions, and to avoid 
any possible interference due to magnetic fields, we have re- 
moved the middle third of the spectral images for the fi » 0.4 




Fig. 2. Equivalent width vs. ju for the oxygen Unes, for different observations. Last panel shows the three 777 nm Oi triplet lines 
(from strongest to weakest: 777.19 nm, 777.41 nm, 777.53 nm). The smooth variation of the equivalent widths with yU in our data 
indicates a low statistical error and highlights its high quality. Equivalent widths for Allende Prieto et al. ( 2004] l were computed 
with our analysis routines using their spectra. 



set, so that the final spectrum was averaged from quiet regions 
only. This procedure reduces the S/N for this set by about 20%. 

After obtaining the spatial and temporally averaged spectra 
for each fi the continuum levels are found manually. In some 
cases like the Oi 615.81 nm line, a local continuum is used 
to compensate for blends not included in the line synthesis. 
The wavelength calibration, detailed in Paper I, is linked to the 
Fourier Transform Spectrograph (FTS) disk-centre intensity at- 
las of jBrault & Neckel ( |1987[ ). It consists on identifying the same 
spectral lines in the FTS atlas and our observations, and then 
obtaining the dispersion relation by a polynomial fit to the line 
cores. It lacks precision, and for the analysis of some lines a 
more precise wavelength calibration is necessary. In the analy- 
sis section we detail a correction to the wavelength calibration 
using the atmosphere models and Fe i lines. The reduced spectra 
for disk-centre and the limb can be found in Fig.[T] 

2.2. Comparison with previous studies 

The oxygen lines we observed have been the object of sev- 
eral studies in the past. Most notably, the centre-to-Umb vari- 
ation of these lines has been studied in Miiller et al.| ( |1968| l, 
|Altrock| ( il968 ) and [Allende Prieto et alT] ( |2004| l. Thelast work 
covers only the 777 nm triplet lines. These works use differ- 
ent instrumentation and techniques, which makes them fairly 
independent of our observations. We compare the equivalent 
widths of the present work with the ones from these studies in 
Fig. |2] Equivalent width measurements have a somewhat sub- 
jective component due to different normalizations and different 
wavelength ranges where it is computed. However, it is the only 
comparison that can be made at least with Miiller et al.| ( |1968| l 
and Altrock ( 1968 ). With [Allende Prieto et al.| ( |2004| i, we used 
the available online spectra and computed the equivalent widths 
in the same way as for our data (direct integration, same wave- 
length range). It should be noted that their normalization around 
the 777.41 and 777.53 nm lines is slightly different from ours, 
leading to a higher continuum level. This in turn makes the lines 
stronger when computing the equivalent widths by direct inte- 
gration. The equivalent widths as a function of fi for all the lines 
included in this work are listed in Table[T] The uncertainties arise 
mostly from the continuum placement. For the lines in the 777 
nm region, the errors are larger because a larger integration re- 
gion was used. 



3. Theoretical line profiles 

We employ the 3D time-dependent hydrodynamical simulation 
of the solar photosphere used in Paper II (Trampedach et al.. 



in preparation; Asplund et al., in preparation). This model was 
computed using a more recent version of the Stein & Nordlund 
codes. When compared with the 3D model of [Asplund et al. 
(2000), it includes more detailed radiative transfer: 12 opacity 
bins, additional sources of continuum opacities such as photo- 
ionization cross-sections from the Opacity and Iron Projects 
( [Cunto et al.|1993|[Hummer et al.|1993] l, etc. In Fig. [3] we show 
the differences in the mean temperature structure (averaged over 
surfaces of constant optical depth) between both 3D models. 

A comparison is also made with one-dimensional, hydro- 
static, time independent solar model atmospheres, including the 
semi-empirical model of Holweger & Miiller ( 1974| l and the 
MARcs model ( [Gustafsson et al.|2008| l. 

We computed the synthetic LTE and NLTE (oxygen only) 
line profiles, making use of our LTE line formation code and the 
multi3d code ( Botnen 1997[ [Asplund et al.[[200 3]l. For a more 
detailed description of the 3D model, line formation and atomic 
data used we refer to Paper I However, a few things are done 
differently from [Paper I[ as discussed below. 

In this work we make use of 90 snapshots of the 3D simula- 
tion (in Paper I we use 20 snapshots). These cover x 45 min of 
solar time and have (red) = 5778 + 14 K. To save computational 
time, 3D NLTE line formation is carried only on four snap- 
shots (chosen to have a representative range of effective tem- 
perature and time coverage). Using multi3d we compute NLTE 
and LTE line profiles, which are used to obtain the wavelength- 
dependent NLTE/LTE ratio for each fi value. In contrast with 
Paper I the NLTE/LTE ratio is averaged over all the spatial 



points and the four snapshots. For each fi value, we then mul- 
tiply the NLTE/LTE ratio by the spatially and temporally aver- 
aged LTE line profiles (computed with our LTE code for all the 
90 snapshots). The NLTE line profiles were computed for eight 
different values ofSn- 0.01,0.1,0.3, 0.5, 1, 1.5, 3, 10. 

The line profiles from the ID models, based on the marcs and 
Holweger-Miiller model atmospheres, were computed using the 
same procedures, opacities and radiative transfer. The exception 
is that a microturbulence of ^tuib - 1-0 kms"' has been used 
when computing these line profiles. In addition, ID line profiles 
were convolved with a Gaussian to account for macroturbulence. 
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Table 1. Equivalent widths for the lines included in this analysis, as a function of /i. 
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W [pm] 










Species 


A [nm] 


yu = 0.999 ± 


0.001 


= 0.793 ± 


0.012 


yci = 0.608 ±0.018 


H = 0.424 ± 


0.024 


fi = 0.197 ± 


0.003 


Fei 


615.1618 


4.95 


±0.04 


5.14 


±0.04 


5.19 ±0.04 


5.55 


±0.04 


5.65 


±0.04 


Si I 


615.5693 


0.61 


±0.03 


0.62 


±0.03 


0.68 ± 0.03 


0.67 


±0.03 


0.62 


±0.03 


Cai 


615.6023 


0.92 


±0.03 


0.92 


±0.03 


1.01 ±0.08 


0.97 


±0.08 


1.14 


±0.03 


Fei 


615.7728 


6.45 


±0.04 


6.49 


±0.04 


6.39 ± 0.05 


6.58 


±0.16 


6.42 


±0.04 


Oi 


615.8186 


0.63 


±0.03 


0.50 


±0.03 


0.38 ±0.06 


0.33 


±0.03 


0.19 


±0.03 


Fe I 


615.9378 


1.21 


±0.04 


1.28 


±0.07 


1.31 ±0.06 


1.39 


±0.11 


1.39 


±0.04 


Fe I 


629.0965 


7.06 


± 0.05 


7.05 


± 0.05 


6.92 ± 0.05 


6.87 


± 0.15 


6.36 


± 0.05 


Fei 


629.7793 


7.54 


±0.04 


7.53 


±0.04 


7.52 ±0.15 


7.77 


±0.04 


7.71 


±0.04 


[Oi] 


630.0304 


0.45 


±0.04 


0.50 


±0.04 


0.54 ± 0.04 


0.58 


±0.04 


0.65 


±0.04 


Sen 


630.0698 


0.59 


±0.03 


0.65 


±0.03 


0.68 ± 0.03 


0.76 


±0.03 


0.77 


±0.03 


Oi 


777.1941 


8.48 


±0.09 


8.00 


±0.14 


7.21 ±0.09 


6.38 


±0.10 


5.09 


±0.10 


Oi 


777.4166 


7.44 


±0.09 


6.96 


±0.15 


6.21 ±0.09 


5.44 


±0.10 


4.22 


±0.33 


Oi 


777.5390 


5.90 


±0.09 


5.47 


±0.11 


4.83 ± 0.09 


4.17 


±0.10 


3.08 


±0.28 


Fei 


778.0557 


14.81 


±0.14 


14.45 


±0.14 


14.13 ±0.14 


13.39 


±0.14 


11.66 


±0.14 



Mean 3D Model (new) 
Mean 3D Model (old) 




Table 2. Lines studied in the present work and their parameters. 



Fig. 3. Comparison of the mean temperature structure of the new 
3D model used in this work and the 'old' 3D model of |Asplund| 
et al. ( 2000| i, plotted against the optical depth at 500 nm. 



Atomic Species 


A [nm] 
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Fe I 
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-3.299"= 


2.18*= 
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0.92 


Fei 


615.7728'' 
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0.63 
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615.9378'' 
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Fe I 


629.0965" 


-0.774" 


4.73" 


7.06 


Fei 


629.7792"= 


-2.740'= 


2.22*= 


7.54 


[Oi] 


630.0304"= 


-9.717"= 


0.00"= 


0.45 


Sen 


630.0698'' 


-1.887'' 


1.51" 


0.59 


Oi 


777.1944"= 


0.369'= 


9.15*= 


8.48 


Oi 


777.4166"= 


0.223'= 


9.15'= 


7.44 


Oi 


777.5390"= 


0.002'= 


9.15"= 


5.90 


Fei 


778.0557'' 


0.029'' 


4.47" 


14.81 



Equivalent widt hs at disk-centre measured from our obser vations 
VALD database jPiskunov e t al.'l995l, |Kupka et al.|2000 1 
NIST database dRalchenko et al. 2009) 



Naveetal. dl994i 



Storey & Zeippen| ( [2000l i 



The adopted macroturbulences vary from line to line, as detailed 
in Sect.|4] For consistency, ID NLTE line formation is computed 
in a similar manner to the 3D model, obtaining the NLTE to LTE 
ratio from multi3d and multiplying it by the profiles obtained 
with the LTE code (even though for the ID models there is only 
one snapshot). 

Our main focus are oxygen lines. However, a few lines from 
other elements are also present in our observations. These lines 
and their properties, along with the oxygen lines, are listed in 
Table [2] In Sect. 4.4 we briefly discuss the results for the centre- 
to-limb variation of these lines. 

For the oxygen lines and, when available, for lines of other 
elements the collisional (van der Waals) broadening was com- 
puted using the quantum mechanical theory of Anstee & O'Mara' 
( [1995 ); Barklem & O' Mara ( 1997 ); Barklem et al. ( 1998 ), avoid- 
ing the need for conventional Unsold ( 1955 1 enhancement fac- 
tors. 



4. Results 

The synthetic line profiles are compared with the observations by 
means of fitting and equivalent width. Before a comparison was 
made with our observations, the synthetic line profiles were con- 
volved with a Gaussian (equivalent to A/AA-2QQ 000) to account 
for the instrumental profile of the spectrograph. When compar- 
ing line profiles with the FTS intensity atlas no convolution took 
place, due to the very high resolution of the atlas. The line pro- 
file fitting was performed using the L-BFGS-B optimization al- 
gorithm ( |Byrd et al.|l995[ l 



4.1. 01 infrared triplet 
4.1.1. Context 

The three O i lines around 777 nm are strong and lie in a part 
of the solar spectrum relatively free of blends and telluric lines. 
These two characteristics alone make them a good abundance 
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O I 777.19 nm O I 777.41 nm O I 777.53 nm 




M M M 

Fig. 4. Top panels: equivalent width vs. ju for the three triplet lines, using 3D LTE and NLTE with different Su coefficients for 
hydrogen collisions. The oxygen abundance was adjusted so that the models have the same equivalent width as the observations at 
disk-centre. Bottom panels: difference in fitted abundance from disk-centre to a given position in jj.. Synthetic profiles were fitted 
against the observations at each /i, and the figures show how much the fitted abundance varies from the fitted value at disk-centre. 



indicator However, these lines show significant departures from 
LTE in the Sun ( |Altrock||1968| ), because of a radiation field 



O I 777.41 nm MARCS 



O I 777.41 nm Holweger-Miiller 



weaker than Planckian in the line formation region (e.g. Eriksson 
|& Toft|1979[|Kiselman|1993l ). 



Aside from the increased computational expense of comput- 
ing the 3D NLTE radiative transfer for these lines, departures 
from LTE introduce additional uncertainties stemming from the 
input physics, in particular photo-ionization rates and coUisional 
cross-sections with electrons and Hi. The latter, as noted in 
Sect. [T] are often taken from the classical estimates and scaled 
by an 5 H factor For the Oi 777 nm lines a scaling factor for 
the Hi collisions was tried by Kiselman ( 1993| l, to reconcile 
different ID models with observed centre-to-limb variation of 
the equivalent widths. Further investigation on 3D and NLTE 



effects on these lines was provided by Kiselman & Nordlund 
( |1995 ), using an early 3D model and showing the feasibility of 
the fine centre-to-limb variations to probe for the atomic param- 
eters (i.e., Sh)- This suggestion was followed by |Allende Prieto| 
|et al.| ( [2004l hereafter AAF04), who used the centre-to-limb vari- 
ation of the O I 777 nm lines and a 3D model to empirically de- 
duce that 5h = 1 was preferrable to 5h = 0. 
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Fig. 5. Same as Fig.|4]but for the ID models and only for the O i 
777.41 nm line. 



Our work differs from lAAF04l in that we use a different set 
of observations and 3D model, a greater range of 5h values, a 
more up-to-date atom and more detailed radiative transfer for 
the NLTE calculations (more simulation snapshots included in 
the full 3D NLTE, more angles used in multi3d). 



4.L2. Comparison with observations 

Our goal is to find the Sw that best describes the centre-to-limb 
variation of strength of the lines. For this purpose we have two 
diagnostics: equivalent widths and line profile fitting. 
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Measuring the variation of line strengths by profile fitting has 
the advantage of being less sensitive to systematics (e.g. blends, 
noise) than the equivalent widths. But when the shape of the 
line profiles does not match the observations, the line profiles 
adjusted in abundance will have a slightly different equivalent 
width than the observations. We present the results using both 
diagnostics, and they both indicate a very similar result. 

To fit the line profiles we vary two parameters: oxygen abun- 
dance and wavelength shift (of the observations). The latter is 
necessary because of the uncertainty in the wavelength calibra- 
tion. For the ID models another free parameter was allowed in 
the fit: macroturbulence. Ideally it should not be allowed to vary 
freely, but rather extracted from nearby lines. The scarcity of 
lines in our observed window around 777 nm makes this task 
difficult; in our observed window there is only one other line, a 
strong Fe i line. Deriving the macroturbulence from only one line 
would probably introduce a similar or larger error as allowing it 
to be a free parameter in the fit for the 1 777 nm lines. 

The results for profile fitting and equivalent width for the 3D 
model are given in Fig. |4] For the equivalent widths, the oxy- 
gen abundance was adjusted for each line so that it matched 
the observations at disk-centre. The same abundance was then 
used for all the other values of jj. Both diagnostics indicate that 
S'h = 1 gives the best agreement with the observations, with a 
small variation from line to line. It is shown with a high confi- 
dence level that LTE line formation is not a valid approximation 
for these lines. In Table|3]we list the derived abundances from fit- 
ting the disk-centre line profiles. They were fitted in the follow- 
ing wavelength ranges: 777.175-777.24 nm, 777.36-777.46 nm 
and 777.50-777.565 nm. 

Corresponding results for the ID models are given in Fig.|5] 
but only for the Oi 777.41 nm line; the other lines behave 
similarly. They indicate a mixed scenario. On one hand, the 
Holweger-Miiller model gives a similar result to the 3D model: 
the best agreement is with 5h = 1. But for the marcs model no 
value of 5h can give a reasonable agreement with the observa- 
tions, both in equivalent width and line profile fitting. 

In Fig.|6]the 3D model predicted line profiles for disk-centre 
(adjusted in abundance) and the limb (using the disk-centre 
abundance), are shown for LTE and Sh = 0.01, 0.3, and 1. It can 
be seen that LTE performs very poorly at the limb (too weak), 
that Sh = 0.01 does not fare much better (too strong) and that 
5h = 1 shows a much better agreement. 

One can also see that at disk-centre the LTE profiles seem to 
fit the observations better The disk-centre profiles for 5h = 1 are 
narrower and deeper than the observed. At the limb, in terms of 
shape only, = 1 profiles have a much better fit. The Sh that 
best describes the centre-to-limb variation of the line strengths 
does not seem to be the best at describing the shapes of line 
profiles at disk-centre. This same effect is more obvious when 
looking at the granulation variation of the FWHM at disk-centre 
Paper I| . It may very well be connected with the finding of 
'aper I that no single Sh can provide the best agreement with 
the observed equivalent widths for both granular and intergran- 
ular regions. The reason for these discrepancies is yet not clear. 

To better quantify the agreement with the observations of 
different Sh, ^ minimization was made with the equiva- 
lent widths, simultaneously for the three lines and for the five 
jU values. For each value of Sh we varied the oxygen abun- 
dance and found the value that minimizes the squared difference 
between observed and predicted equivalent widths, weigthed 
by the observational error bars. The reduced is defined as 
i/N ■ Yj (Wobs — VKmodei)^ /o"^, whcrc is the number of degrees 
of freedom (in this case, - 15-1). Thex^ values, as a function 



Table 3. Derived oxygen abundances (log eo) from fitting disk- 
centre line profiles from the 3D model, using 3D NLTE radiative 
transfer for different values of Sh both for our SST observations 
and the FTS atlas of |Brault & NeckeT| ( fr987] l. 



1 line [nm] 


Sh 


SST/TRIPPEL 


FTS atlas 


777.19 


0.01 


8.50 


8.52 




0.1 


8.53 


8.55 




0.3 


8.58 


8.60 




0.5 


8.61 


8.63 




1 


8.66 


8.68 




1.5 


8.73 


8.75 




3 


8.76 


8.77 




10 


8.82 


8.84 




LTE 


8.87 


8.88 


777.41 


0.01 


8.54 


8.55 




0.1 


8.56 


8.57 




0.3 


8.60 


8.61 




0.5 


8.63 


8.64 




1 


8.69 


8.69 




1.5 


8.74 


8.75 




3 


8.76 


8.77 




10 


8.82 


8.83 




LTE 


8.85 


8.86 


777.53 


0.01 


8.58 


8.58 




0.1 


8.60 


8.60 




0.3 


8.63 


8.64 




0.5 


8.65 


8.66 




1 


8.69 


8.70 




1.5 


8.74 


8.75 




3 


8.75 


8.76 




10 


8.80 


8.81 




LTE 


8.83 


8.83 



of 5h are shown in Fig.|7] Comparing for all the 5h we identify 
the value that gives the best agreement, using a parabolic fit to 
the lower values. 

For the 3D model we find the minimum to be at 
Sh ~ 0.85. We use this value to derive the oxygen abundance 
from these lines. Averaging over the three lines we obtain an 
oxygen abundance of log eo - 8.68 for the 3D model. Repeating 
the procedure for the ID Holweger-Miiller model we also obtain 
Sh ~ 0.85, deriving an abundance of 8.66. Although with the ID 
MARCS model no 5h reproduces the observations, for comparison 
we derive an abundance of 8.61 using Sh = 0.85. 



4.L3. Comparison with previous worl< 

At first glance our findings are consistent with the Sh - 1 es- 
timated by AAF04I But these two findings cannot be directly 
compared, as we use a different 3D model. To compare our 
methods and observations with |AAF04| we repeated our anal- 
ysis for the same 3D model they use. The best fitting value is 
then Sh ~ 0.3, which gives a similar agreement with the obser- 
vations as S'h = 1 with the new 3D model. Between Sh - 
and S'h = 1 (the two values tested by AAF04| l we find that the 
latter gives a better agreement with the observations, consistent 
with the findings of|AAF04[ But the best fit of Sh ~ 0.3 indi- 



cates a model dependence in this derivation, also seen for the 
ID models. Interestingly, although the best fitting Sh varies for 
the two 3D models, the abundances derived from profile fits with 
the best choice 5 h are almost unchanged between the two mod- 
els (^ 0.01 dex). 
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Fig. 6. Our observations (dots) and 3D synthetic line profiles (solid line). Each panel shows the 777 nm triplet lines at disk-centre 
and // - 0.2, for LTE and three different values. The line profiles at disk-centre were fitted to the observations, while at fi - 0.2 
the abundance was set fixed to the disk-centre abundance, while the Waveshift was adjusted. For clarity, disk-centre line profiles 
have been subtracted 0.05. 



4.1.4. Effect of blends and electron collisions on the inferred 
Sh and abundances 

There are a few weak CN and C2 Unes in the 777 nm region. The 
VALD database also lists some weak atomic lines. Most of these 
lines are very weak and their effect on the O i lines is negligible. 



At disk-centre the known blendaHof the Oi lines are very faint. 



' There are a few unknown blends visible in the disk-centre spectrum. 
IVIost notably, the feature at the 'knee' of the wing of the 777. 19 nm line. 
Their effects in the fit were minimized by limiting the profile fitting 
region to parts not visibily affected by blends. 
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But as molecular lines get stronger at the higher layers probed 
by the limb spectra, it is important to determine if they have an 
effect on the inferred (the effect on disk-centre abundance 
obtained by profile fitting is negligible). To investigate for such 
effects we have computed the spectra of 1 777 nm over a larger 
spectral region and including additional molecular and atomic 
lines, with data from J. Sauval (2008, private communication). 
Our testing showed that the blends have very little influence in 
the derived Sh value. Because the atomic and molecular data for 
these weak blending lines is not known accurately, they have not 
been included in our main analysis. 

We also investigate the effect of computing the O i lines si- 
multaneously: their wings are very extended and overlap, es- 
pecially in LTE and for 777.41 and 777.53 nm. For LTE at 
disk-centre, computing the three O i lines simultaneously caused 
a reduction of 0.002 dex in the fitted oxygen abundance of 
777.41 nm line. At the limb there is virtually no difference be- 
tween computing the lines separately or together, because the 
lines are much weaker For the same line strength the LTE pro- 
files have the most extended wings when compared with NLTE 
(the extent of the wings gets smaller with smaller S h) and thus 
it is for the LTE profile that the issue of overlapping wings gets 
more important. Our single line approximation is thus fully jus- 
tified. 

In this work we make use of revised electron collision data 
( [Barklem 2007) . Although this data is believed to be the most ac- 
curate currently available, it is important to check to what point 
our empirically estimated depends on the adopted electron 
collision rates. To assess this effect, we have performed several 
tests by artificially changing the electron collisions. Due to time 
constraints these tests have been performed in ID models, al- 
though they were confirmed with a few 3D runs. An interesting 
find is that if the electron collision rates are artificially decreased, 
the best fitting S h value is still the same (even if we neglect elec- 
tron collisions altogether). Obviously the fitted abundance values 
will differ (especially for Sh < 1), but the best fitting value re- 
mains the same. This result indicates that H i collisions are more 
important than electron collisions in the centre-to-limb varia- 
tion of the lines, because the electron pressure decreases more 
quickly with height. On the other hand, if the electron collision 
rates are increased by one order of magnitude it will change the 
best fitting Sh value. This makes sense, as an increase on the 
collisions from electrons will shift the result more towards LTE, 
regardless of the 5 h used. 

4.2. [O I] 630.03 nm and Ni I blend 
4.2.1. Context 

Although the forbidden O i line at 630.03 nm is the strongest of 
the [Oi] lines it is still very weak in the Sun. Unlike the per- 
mitted lines, LTE is a very good approximation in the formation 
of the [O i] line (see Asplund et al. 2004). Accordingly, we only 
use the synthetic LTE profiles in the following analysis. Being a 
weak line it is highly susceptible to even faint blends. In partic- 
ular, it is blended with a Ni i fe ature at 630.0335 nm (|Lambert| 
[T978)[Ane nde Prieto et al.|2001| [Johansson et al.|2003i For the 
Sun, the contribution of the Ni i component to the blend is about 
30%, depending on the Ni abundance used. With such a signifi- 
cant blend it is of paramount importance to properly account for 
the Ni contribution. This work is similar to the work of' Allendel 




Fig. 7. Reduced as a function of Sh from the minimization of 
the equivalent widths (see text), for the 3D model. The LTE value 
is indicated (dashed line). The minimum at 5h = 1 indicates the 
best agreement with the observations. Inset: same plot in more 
detail for 5 h < 2, showing the parabolic fit to the lower values 
(thick solid line). 

Table 4. Ni i isotopes used and scaled log gf values. Wavelength 
of main component and isotopic splitting from Johansson et al. 
(2003) , isotopic fractions from Rosman & Taylor, ( ,1998) . 



Ni isotope 


% 


Wavelength [nm] 




58Ni 


68.27 


630.0335 


-2.2757 


«ONi 


26.10 


630.0355 


-2.6933 


61 Ni 


1.13 


630.0365 


-4.0569 


62Ni 


3.59 


630.0375 


-3.5549 




0.91 


630.0395 


-4.1509 



also'Caf fau et al.| (|2008) and 'Ayres' (2008) have performed simi- 
lar studies using a different 3D model. 

4.2.2. Nil blend 

Properly accounting for the Ni i blend means one needs to know 
its properties with high precision. Johansson et al. ( 2003| l have 
measured the log gf for this Ni i transition to be -2. 1 1 . Because 
of the presence of several Ni isotopes and their isotopic splitting, 
the Ni contribution at 630.03 nm is given as a series of blends. 
In our treatment we included the five most abundant Ni isotopes 
scaling the log by their isotopic abundance (see Table |4]l. We 
have used a fixed Ni abundance of log eNi - 6.22, as inferred 



from other Ni lines using the present 3D model (Asplund et al. 
|2009l ). 

An important difference between the present work and 
lALAOll was the treatment of the Ni i blend. We used a fixed 
Ni abundance of 6.22 and adopt the measured \og gf of -2.11, 
whereas at the time of |ALA01| there was no experimental 
\og gf and Ni was treated as a free parameter. Their best fit- 
ting log(g/eNi) - 3.94 is equivalent to a Ni abundance of 6.05 
with the present value of log gf, smaller than the Ni abundance 
we use. In ALAOl the blend of the Oi and Nil was done by 



Prieto et al. (2001 hereafter ALAOl), that included an analysis 
of [O i] 630.03 nm with a 3D model (albeit an earlier version), 
but for the flux profile instead of intensity at different fi values; 



co-adding the individual flux profiles for the different species. 
Because these are very weak lines, that was a reasonable approx- 
imation. However, the approximation of co-adding the profiles 
starts to break down as one looks at profiles of smaller /i values, 
where the difference between co-adding and a proper treatment 
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of the blend (summing the opacities for each component) can 
amount to a few percent of the equivalent width. In this work 
we opted to treat the blend properly, instead of co-adding the 
profiles. 

4.2.3. Wavelength calibration 

Another important point in the treatment of this line is the wave- 
length calibration. As Ayres (2008) noted, the balance between 
the components of the blend is sensitive to velocity errors of a 
few hundred ms or a few mA. If the wavelength calibration 
translates the line to be blueshifted from its 'true' location the 
fitting will favour a higher oxygen abundance. There is a degen- 
eracy between the wavelength shift and the Ni abundance, which 
makes it difficult to determine the best-fitting values when both 
are allowed to vary. 

Unfortunately the absolute wavelength caUbration of neither 
our observations nor the FTS atlas is precise enough for the anal- 
ysis of this blend. Hence we adopt a correction to the labora- 
tory wavelength similar to that of Ayres| ( 2008[ l: using strong 
Fei lines in the neighbourhood to find the wavelength differ- 
ence between the 3D model and the observations. By measur- 
ing the wavelength difference between the synthetic profiles and 
the observations for a set of reference Fe i lines, one can effec- 
tively put the profiles in the same 'laboratory' wavelength scale. 
This approach has the advantange of removing the systemat- 
ics associated with the observations {e.g. solar rotation) and any 
shortcomings in the 3D model (e.g. overshooting of convective 
blueshifts). We chose the three Fei lines at 629.77, 630.15 and 
630.24 nm - as they were available in the limited wavelength 
window of our observations - to calibrate the wavelength scale. 
To measure the wavelength difference we fit the theoretical ID 
and 3D line profiles of these Fe i lines and allow a wavelength 
shift of the observations as a free parameter. The mean wave- 
length shift from the three lines was extracted for line profiles at 
all the solar disk positions of our observations. The standard de- 
viation between the shifts from the three Fe i lines were 19, 27, 
33, 51, 140ms-i for respectively = 1,0.8,0.6,0.4,0.2. Except 
the yu = 0.2 measurement, the numbers are an encouraging sign 
of how precise the wavelength calibration is. For the observa- 
tions at fi = 0.2, due to the scatter between Fe i line indicators 
we adopted the wavelength shift as given by the Fe i 629.77 nm 
line only. 

If the wavelength shift of the synthetic line profile was al- 
lowed to be a free parameter in the fit, then the best fitting shift 
to the FTS intensity atlas is very close (< 75 m s"') to the fixed 
shift obtained by the Fei lines procedure. Its effect on the de- 
rived abundance is thus small (< 0.01 dex). Even if one accounts 
for some uncertainty around the wavelength shift, the effects for 
the oxygen abundance are small. A relatively large variarion of 
(+330 m s"') around the used Waveshift for the FTS intensity 
atlas has an effect on the oxygen abundance of x +0.015 dex, al- 
thought the agreement with the observed profile deteriorates sig- 
nificantly. These variations of the oxygen abundance are much 
smaller than the ones found by |ALA01 when trying different 
shifts. The difference is likely to be due to the fact that we kept 
the nickel abundance fixed, while | ALAO 1 1 allowed it to be a free 
parameter in the fit. With nickel fixed, the fitting procedure can- 
not compensate with more or less nickel, which means that the 
line shifting has a smaller effect on the oxygen abundance. 

Besides the absolute wavelength calibration of the observa- 
tions, there is also the issue of the uncertainties in the laboratory 
wavelengths for the oxygen and nickel transitions. The total un- 



certainty in the wavelength difference between the lines is ap- 
prox. ±0.236 pm, being ±0.2 pm for [O i] ( Eriksson|1965 \ and 
±0. 125 pm for Ni i (Joha nsson et al.|2003 i, weighted by the iso- 
topes. We tested the Tine profile fitting for the lines using the 
two extreme cases of wavelength differences allowed by the er- 
rors and found that its effect on the derived oxygen abundance 
is negligible (< 0.01 dex). Once again, this value is smaller than 
the 0.05 dex found by ALAO 1 , because we do not have the nickel 
abundance as a free parameter. 



4.2.4. Continuum level and fitting range 

The local continuum of the [Oi] + Nil line wings is slightly 
depressed by the presence of nearby Fe i and Si i lines, and also 
about 70 weak CN lines between 630.00-630.06 nm, as noted 
by |ALA01| To account for this effect, we multiplied each set 
of observations by a factor C. For the FTS intensity atlas C = 
1 .0054 was used. It is important to note that the determination of 
this local continuum is a major source of the total uncertainty. If, 
for example, one just took the FTS intensity atlas's continuum 
and computed the equivalent width for this line, the equivalent 
width would be a: 30% higher. Care was taken in determining 
the local continuum for our observations, but one should keep in 
mind the uncertainty attached to it. 

Another important but seldom discussed point in the profile 
fitting is the wavelength range where to fit the profile. 'ALAOl 



have done the profile fitting in the range 630.015-630.04 nm, 
others are not so explicit. Ideally one would like to perform the 
profile fitting in a range as broad as possible to include the max- 
imum amount of information about the line. The presence of 
blends, especially unknown features, forces one to narrow the 
fitting range to avoid the influence in the fit of features not in- 
cluded in the line synthesis. We have adopted the same wave- 
length range as ALAOl^ as it avoids several weak blends but 
includes most of the [Oi] line. 



4.2.5. Comparison with observations 

The line profiles for the [O i] -H Ni i blend were computed for 
a fixed log fNi and assuming LTE. In the fit for this line pro- 
file we only allowed one free parameter: the oxygen abundance. 
The Waveshift for each position was obtained by the Fe i lines 
method, and the local continuum multiplication factor C ob- 
tained empirically, as noted above. 

We present the results for the equivalent width analysis over 
the observed positions in the solar disk, the profile fitting of our 
observations at /u - 1, 0.2 and the FTS intensity atlas for com- 
parison with other studies. Some results from line fitting are in 
Fig. [8] and the equivalent width results are in Fig. |9] Fitting the 
line profile to the FTS intensity atlas we obtain an oxygen abun- 
dance of log eo - 8.66, while for our observations at jj. - I we 
obtain 8.64. In the variation of equivalent width with fj. in Fig.|9] 
one can see that the 3D model predicts a slightly stronger line as 
fi decreases. In terms of fitting the line profiles for different yU the 
variation in fitted abundance for the 3D model is never more than 
0.02 dex in regards to the disk-centre abundance. In the right- 
most column of Fig.|8]we plot the observations atjj- 0.2 with a 
synthetic profile using the oxygen abundances derived from the 
disk-centre fits (our observations), the continuum level adjusted 
independently. The 3D model at /i = 0.2 is slightly stronger than 
the observations, and the line profile is slightly wider and not as 
deep as the observations. Still, in terms of best-fitting abundance 
the difference is minimal. 
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3D Model: FTS Intensity Atlas 
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Fig. 8. Synthetic line profiles (solid lines) for the O i+Ni 1 630.03 nm blend vs. observations (circles). Shaded regions indicate range 
over which the profiles were fitted. Left panels: fit for FTS Intensity atlas, done over the wavelength range of Allende-Prieto et al. 
(2001). Middle panels: fit for our observations, disk-centre. Right panels: our observations at = 0.197, and profiles from models 
using abundance fitted at disk-centre, adjusted for continuum. The wavelength shift was measured for each model and disk position 
separately from Fei lines (see text). Ni abundance was fixed at 6.22, 6.16, 6.26 for the 3D model, ID MARCS model and ID HM 
model, respectively. 



For the ID models we have applied the same wavelength cor- 
rection with Fe I lines, repeating the procedure for line profiles 
synthesized with the ID Holweger-Miiller and marcs models. 
The best fitting macroturbulence was also extracted from these 
Fei lines, and then used for Oi+Nii. The Holweger-Miiller 
(HM) model does a good job at describing the line strengths, as 
evidenced by the equivalent widths in Fig.|9] However, a closer 
inspection at the line profiles shows a red wing that is too strong. 
This may possibly indicate an excess of nickel: for this model 
log fNi = 6.26 was used, as derived with the same model using 
other Ni lines. Looking at Fig. [8] in particular at the residuals, it 
is clear that the agreement is not as good as for the 3D model. 
This is expected, as the ID models do not reproduce the convec- 
tive motions in the photosphere as well as the 3D model, even 
using micro and macro turbulences. In particular for the HM at 



fi - 0.2 the macroturbulence derived from the Fe i lines is prob- 
ably low, as the line profile is noticeably deeper than the ob- 
servations, while the equivalent widths are very similar. For the 
MARCs model the trend in equivalent width vs. /i is different: the 
model predicts slightly weaker line strengths than the observa- 
tions. The fitted line profiles at disk-centre also indicate an ex- 
cess of intensity in the red wing. The line profile at low /u with 
the disk-centre abundance shows a slightly weaker line than the 
observations. The oxygen abundances from profile fitting with 
the FTS disk-centre atlas for the HM and imarcs models were, 
respectively, 8.69 and 8.66. 
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Fig. 9. Equivalent width vs. fi for the [O i]+Ni i 630.03 nm lines 
blend, for our observations and different models at different po- 
sitions in the solar disk. 



4.2.6. Comparison with previous work 



An easy comparison to be made is with |ALA01| Some things 
have been done differently in the present study. We compute a 
proper blend of the two transitions (adding line opacities), in- 
stead of computing separate line profiles for each line and co- 
adding them. We also perform a more precise wavelength cali- 
bration usinging Fe i lines. And finally, we employ a more recent 
3D Model and a higher Ni abundance. Fitting the flux profiles 
against the FTS flux atlas and using the same fitting range in 
ALAOl we obtain an oxygen abundance 0.07 dex lower. The 



reason for this difference is mainly due to a higher Ni abun- 
dance used (6.22 vs. 6.05), which translates to a s; -0.07 dex 
difference in the derived oxygen abundance. The proper blend- 
ing of the line instead of co-adding profiles causes a difference 
of -0.01 dex in oxygen. The proper wavelength calibration us- 
ing Fe I lines instead of only correcting for gravitational redshift 
causes a difference of -i-O.Ol dex. 

Fitting the line profiles with log eni = 6.22 instead of 6.05 
leads to a better agreement with the observed disk-centre profile 



6.22 



profile Cr^^2 



0.53x1 Q^), ^ worse agreement for the observed flux 



1-514.05)- 



4.3. 01615.81 nm 
4.3.1. Context 

The 1 615.81 nm line is a triplet, with components at 615.8149, 
615.8172 and 615.8187 nm, according to the NIST database. 
The latter is the strongest component. The components are so 
close that they are unresolved in the solar spectrum. Its position 
in the solar spectrum is far from optimal. The line is close to a 
strong (x! 4.8 pm) Fei line at 615.77 nm and there are multiple 
weak blends around. The VALD database lists several atomic 
lines in this region, but besides the strong Fei line only the 
Ndii 615.782 nm. Nil 615.800 nm and Fei 615.803 nm lines 
have any measurable effect, though still just barely noticeable. 
Spectral synthesis using only these lines is still not enough to 
explain the solar spectrum. However, there seem to be several 
molecular lines in this region. These are mostly CN lines, with 
a few C2 lines. The molecular lines are important for our anal- 
ysis and their effect is discussed below. It should be noted that 
even accounting for the molecules there are still unidentified fea- 



tures in this spectral region, most notably the feature at around 
615.793 nm. 

The presence of these blends depresses the continuum level 
around the O i lines, so that it becomes a significant source of un- 
certainty. Our analysis, either in computing the equivalent widths 
from the observations or fitting the line profiles, suffers from the 
same problems. For example, if one assumed this oxygen line to 
have no blends and took the continuum level of the FTS disk- 
centre intensity atlas, then the equivalent width could be up to 
33% higher than our estimate for the oxygen-only contribution. 
This effect is amplified as one looks at spectra closer to the limb 
because the different blends vary differently in strength over the 
solar disk, which makes a precise determination of the contin- 
uum level difficult. Unlike the stronger Oi 777 nm triplet the 
predicted NLTE effect is rather small for the O i 615.8 nm lines: 
Asplund et al. (2004) predict it to be -0.03 dex when neglecting 
H I collisions. 



4.3.2. Wavelengtli calibration 

A similar approach to the wavelength calibration of [Oi] 
630.03 nm was used. Because a precise wavelength calibration 
for this line is not as crucial as [Oi], we derive the wavelength 
calibration using only the nearby Fei 615.77 nm line. It is im- 
portant to make the wavelengths of the O i consistent with this 
Fe I line because they lie on its red wing, and a different wave- 
length difference between them would imply a lower or higher 
blend influence. 

Using the reference laboratory wavelengths from the NIST 
database for the Fei 615.77 nm and Oi 615.81 nm lines and 
calibrating using the Waveshift of Fei lines, the synthetic Oi 
profiles are shifted from the observations. This can be seen in 
the left panels of Fig. [TO] The 3D synthetic profile is shifted to 
the blue in regards to the observations, while for the ID mod- 
els the opposite happens. None of them matches the observa- 
tions perfectly. For the 3D case, the shift to the observations is 
about 0.89km s"', or 1.8 pm. This could mean that there is a sig- 
nificant error in the laboratory wavelengths; the predicted line 
shifts from the 3D models are inconsistent for the Fe i and O i 
lines; or there are blending features causing a shift in the profile. 
Including known molecular and weak atomic lines is not the ex- 
planation. In the top middle panel of Fig. [TO] we have artificially 
shifted the O i lines, so that the synthetic profile matches the FTS 
intensity atlas. For the rest of the analysis we will be using these 
artificially shifted profiles (both for 3D and ID models), but the 
derived abundance from fitting the disk-centre profiles remains 
essentially the same if the profiles are not shifted (but as seen 
in Fig. 



10 the agreement with the observations is worse). The 



results for the equivalent widths are not affected by this shift. 



4.3.3. Molecular blends and comparison with observations 

We have included several CN and C2 lines in the 615.8 nm re- 
gion (J. Sauval, 2008, private communication). The line data are 
given in Table|5] These data were obtained from from J. Sauval 
(2008, private communication). The molecular line data is un- 
certain, especially wavelengths for CN lines. With this uncer- 
tainty in mind, we have adjusted some of the wavelengths and 
log gf of the CN lines from Table |5] The two (9, 4) band lines 
at 615.8522 and 615.8796 nm have been red shifted by 7.2 and 
11.6 pm, respectively, so that they match two features in the so- 
lar disk-centre spectrum. Furthermore, the log gf values of these 
two lines have been increased so that they fit reasonably the ob- 
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3D Model: FTS Intensity Atlas 
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Fig.lO. Synthetic LTE line profiles for the Oi 615.81 nm lines (including blends from Fei, one C2 and five CN lines). All panels 
except lower left show synthetic profiles (solid) vs. observations (circles). Top left panel: location of the Oi line in respect to the 
observations, using Waveshift from Fez lines and laboratory wavelengths from the NIST database. Middle panels: disk-centre line 
profiles for the 3D and ID Holweger-Miiller models with the Oi lines artificially shifted to the red. Right panels: limb line profiles 
for the 3D and ID Holweger-Miiller models, using the same shift as in the middle panels. Disk-centre profiles were fitted for 
abundance, while for the limb profile the disk-centre abundance was used. Lower left panel: synthetic profiles illustrating the effects 
of the atomic and molecular blends at disk-centre and limb. Shaded regions indicate range over which the disk-centre profiles were 
fitted. 



served features. For consistency the logg/ values of the other 
three CN lines have been increased by the same amount. For 
the 3D model this increase in logg/ amounts to 0.35 dex; we 
note that CN lines from different bands do not necessarily re- 
quire the same correction in abundance or log gf . However, we 
used the same log gf correction because the data for the differ- 
ent CN lines have similar uncertainties. Additionally, the (10, 4) 
band line at 615.8223 nm has also been red shifted by 9.3 pm. 
As detailed below, the inclusion of molecules is of paramount 
importance to understand the solar spectrum at the limb. 

Our synthetic profiles included a total of ten blends: the three 
Oi lines, the strong Fei Une and the molecular lines detailed 
above. There are other atomic lines in the region, not included 
because they are very weak. The profile fits in Fi g. [T0| include 
all these blends. In the lower left panel of Fig. 10 is shown 
the synthetic spectrum of only O i and Fe i compared with the 
molecules, at the disk-centre and limb. At 10.741 eV, the oxygen 
lines have a high excitation energy, making them very weak at 
the higher layers of the atmosphere probed by the limb spectra, 
when compared to disk-centre. The molecules, on the other hand 
are formed in cooler regions and are thus significantly stronger 
at the limb than at disk-centre. As is visible in the figure, at disk- 
centre the molecules blended with oxygen have a very small 
contribution to the total line strength, yet at the limb most of 
the line strength comes from molecules. This type of blending 
can explain the observations well. As seen in the top right panel 



Table 5. Molecular lines included in the 615.8 nm region. 



Mol. 


Wavelength 


System 


Band 


logg/ 






[nm] 








[eV] 


CN 


615.7850 


A-X 


(7, 2) 


-1.498 


1.800 


Co 


615.8106 


rf^Hj - a^n„ 


(0, 2) 


-1.652 


0.508 


CN 


615.8109 


A-X 


(6, 1) 


-2.140 


1.093 


CN 


615.8223 


A-X 


(10, 4) 


-1.613 


2.066 


CN 


615.8522 


A-X 


(9, 4) 


-1.959 


1.062 


CN 


615.8796 


A-X 


(9, 4) 


-1.924 


1.088 



of Fig. [To] the line profile shape and strength using the oxygen 
abundance fitted at disk-centre is very close to the observed pro- 
fil^ Because the molecular lines to the blue of the oxygen line 
dominate at the limb, the line centre is significantly blueshifted 
as is also confirmed by the observations (the wavelength cali- 
bration using the strong Fe i line proved important to establish 
this). Furthermore, the variation in equivalent width in Fig. 
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another confirmation that models including the molecules repro- 
duce the observations satisfactory. 

The fitted oxygen abundances with the disk-centre FTS at- 
las were: logeo = 8.62 (3D Model), 8.72 (ID HM), 8.61 (ID 
iviARCs). For our own observations at disk-centre the abundance 



' The oxygen abundance was the only adjusted parameter in the disk- 
centre fit, and all other lines have the same abundance and parameters 
throughout the analysis. 
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Fig. 11. Equivalent width vs. jj. for the three 1 615.81 nm lines, 
including the Fe i, CN and C2 blends. 



values were higher by about 0.05 dex, difference which can be 
attributed to the difficulty in finding the continuum level. The 
carbon and nitrogen abundances used to compute the molecu- 
lar abundances were different for the 3D and ID models, and 
consistent with the abundances derived from each model from 
other molecular features. For ID models the macroturbulence 
was derived from nearby Fei lines, and a microturbulence of 
^tuib - 1.0 km s"' was used. 

The molecular blends can help to explain the centre-to-limb 
variation in this spectral region. As seen in Fig. 1 1 both ID and 
3D models can explain the variations similarly well (although 
the fitted abundances at disk-centre vary). However, the uncer- 
tainty of the molecular line data and the wavelength shift of oxy- 
gen make this agreement less convincing. These uncertainties do 
not allow us to test reliably the line formation across the solar 
disk, and the effects of different models. Like for 1 777 nm, we 
also experimented with NLTE effects and different Sn. But for 
the 615 nm lines the NLTE effects are so small that with all the 
mentioned uncertainties, it is not possible at this time to extract 
any reliable conclusions about for example the H i collisional ef- 
ficiency. 

4.4. Lines from other species than O I 



In Fig. 12 we present the centre-to-limb variation of the equiv- 
alent width for a few non-oxygen lines included in our ob- 
servations. The properties of these lines are listed in Table |2] 
Nearby blends have been included for the Sii 615.57 nm and 
Cai 616.60 nm. Their wings influence each other, so both of 
them have been included in each of their line profiles, as well as 
a nearby strong Sii line at 615.51 nm. For the Cai 616.60 nm 
we also include the very weak Oi triplet at 615.60 nm, which 
influences the line slightly. Data for all the extra blending lines 
were extracted from the VALD database. 

While the agreement between the 3D model and the lines in 
Fig. 12 is generally worse than for the oxygen lines it is in most 



cases an improvement over the ID models. For most of the Fe i 
lines the ID models predict a much weaker line at the limb, while 
the 3D model prediction usually follows the observations. What 
happens with the Sc 11 630.06 nm is very similar with the figure 
for [Oi] 630.03 nm. This is probably because this line of a ion- 



ized rare element is not very sensitive to the local temperature. 
In the case of Sii 615.57 nm, blends are likely to influence the 
behaviour of the line. This line is in a relatively crowded spectral 
region. For Fe 1 615.93 nm from yU x 0.4 to the limb the ID mod- 
els predict a decrease line strength, while the 3D model predicts 
an increase. The observations in this case seem to follow neither 
and maintain a constant line strength. But it should be noted that 
weak blending features on the wings of Fe 1 615.93 nm make the 
continuum level difficult to determine, hence the big uncertain- 
ties in equivalent width. The Fei 629.77 and 615.16 nm lines 
have a similar variation with fi, which is reasonably reproduced 
by the 3D model while the ID models show larger discrepancies 
and the wrong trend with yU. 

At this time we refrain from drawing conclusions about 
NLTE effects in these lines. However, these may be addressed 
in the future if there is enough atomic data for these atoms. 

5. Conclusions 

The study of oxygen lines over a range of different positions in 
the solar disk has proven very fruitful. From departures of LTE to 
constraints in blends and line data, the detailed line profiles for 
different fi values give us a wealth of information not possible 
before. The relevance of these tests is to ascertain if the temper- 
ature structure of the models can reasonably explain the forma- 
tion of these lines and hence be considered reliable to infer the 
much debated solar oxygen abundance. Although for abundance 
analysis the FTS atlases, given their superior resolution and S/N 
ratio, should be the preferred choice, our high-resolution spec- 
tra have proven decisive to identify the shapes and variations of 
lines and blends over a range of formation temperatures. 

Using detailed 3D NLTE radiative transfer we have com- 
puted the level populations and synthetic profiles for oxygen 
lines. The Oi 777 nm triplet is of particular interest. These lines 
suffer from NLTE effects, the size of which is determined by in- 
elastic collisions with electrons and neutral hydrogen. Due to the 
lack of analytical or experimental data regarding the efficiency 
of the neutral hydrogen collisions, one has to treat it almost as 
a free parameter, which the centre-to-limb variation of the lines 
helps getting an empirical estimate for In a similar approach to 
AAF04 we have determined the empirical value for the multi- 



plication factor of the classical recipe for hydrogen collisions, 
Sh- For the 3D model and the ID Holwer-Miiller model the 
best fit is for Sn ~ 0.85. For the 3D model used by |AAF04| we 
find a smaller value. These differences illustrate the model de- 
pendency of the centre-to-limb derivation of S'h- The ID marcs 
model does not even reproduce the observed equivalent widths 
with any value of 5h- 

A more decisive test of the 3D and other models for the O i 
777 nm lines will have to wait until better data for hydrogen 
collisions is available. Nevertheless, finding an 5h that describes 
the observations provides a self-consistent way of determing the 
extent of the NLTE effects and to extract a reliable value for 
the oxygen abundance. As seen in Tab. [3] the oxygen abundance 
varies quickly with the chosen Sn, especially for 0.3 < 5h ^ 3, 
making a self-consistent estimate of the hydrogen collisions of 
paramount importance. The last remaining discrepancy with the 
empirical Sn approach is that the best fitting NLTE profiles at 
disk-centre are slightly narrower and deeper than the observed. 
It is not clear if this discrepancy can be improved on the model 
side or if it is a consequence of using the classical [Drawin| ( fl968| l 
formulae. 

The weak 615.81 nm Oi line also shows departures from 
LTE, albeit at a much smaller scale than the 777 nm counter- 
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Fig. 12. Centre-to-limb variations of line equivalent widths. The ID models were computed with a microturbulence of 
^turb - 1.0 km s ' and a macroturbulence extracted from nearby Fei lines. The abundance of each element has been adjusted so 
that the models fit the observations at disk-centre. The line profiles of Sii 615.57 nm and Cai 615.60 nm include weak nearby 
blends (see text). All the line profiles were computed assuming LTE. 



parts. Unfortunately the weakness of the NLTE effects and the 
presence of blends make it very difficult to estimate S-a from 
the centre-to-limb variation of this line. When looking at the 
line profiles, especially as one gets closer to the limb, we find a 
line blueshift and strength that are inconsistent with the synthetic 
spectra, when only oxygen is considered. With a high excitation 
energy, synthetic models indicate that at the limb the oxygen line 
should be much weaker than observed. However, if one includes 
a few molecular blends present in the region, the line strength 
and blue shift eff'ects can be almost perfectly reproduced in the 
synthetic spectra. At the limb the molecules blended with the 
oxygen are stronger than the oxygen contribution. Uncertainties 
in the input data for the molecular lines and continuum place- 
ment errors in the observations mean that the precision of the 
results is not enough to discern between the different 5h for the 
weak NLTE eff'ects. 



Blends also play a determinant role in the [O i] 630.03 nm 
line. The unresolvable blend with Ni in this line has been a 
major source of uncertainty in an otherwise good indicator for 
the oxygen abundance (formed in LTE, relatively independent 
of temperature). A good knowledge of the Nil line properties is 
essential. For that we use the recent measurement of its log gf 
value ( [Johansson et al.|[2003| ) and the recent determination of 
the nickel abundance for the 3D model (Asplund et al. 2009 1. 
Establishing the correct contribution of nickel in the blend is not 
enough to guarantee a correct determination of the oxygen abun- 
dance. Because the oxygen and nickel lines are so close together, 
a shift in thek wavelengths will aff'ect the inferred oxygen abun- 
dance, because the fitting procedure will compensate the quan- 
tity of oxygen that best matches the observations. It should be 
noted that because the nickel abundance is kept fixed, the effect 
of the wavelength shift on the oxygen abundance is not as signif- 
icant as if it was allowed to vary {e.g., ALA01| l: we estimate it to 
be less than 0.02 dex. In any case, for this line allowing a wave- 
length shift in the fitting is not desirable, and a precise absolute 
wavelength calibration must be obtained. | Ayres|(200 8 ) made use 
of Fe I lines to calibrate the wavelength scale, as the precision of 
the FTS atlas by itself is not suflicient. In this work we have used 



a similar approach of calibration using Fe i lines to ensure the ac- 
curacy of the results. The line profile fits for the FTS intensity 
atlas using the 3D model show a remarkable agreement between 
model and observations, with less than 0.2% diff'erence for most 
of the line profile. Being not very sensitive to temperature, this 
line displays a small centre-to-limb variation in strength. Even 
with the observational uncertainties inherent of such a weak line, 
we can see that the 3D model provides a very good description of 
the observations. This agreement gives us confidence in the 3D 
model, confirming its robustness to infer the oxygen abundance 
from this line. 

Overall, we find that the 3D model reproduces the observa- 
tions very well. Despite its one-dimensional nature not allow- 
ing for the most accurate shape of the observed profiles, the 
Holweger-Miiller model also reproduces the observed equiva- 
lent widths, which is not surprising because it is semi-empirical. 
The ID theoretical marcs model falls short of its description of 
centre-to-limb variation of the Oi 777 nm and [Oi] 630.03 nm 
lines, probably because its temperature structure departs from 
the solar temperature structure. 

Using the best fitting Sa- 0.85 and abundance results from 
disk-centre line profile fits we infer an oxygen abundance of 8.68 
for the 3D model using the 777 nm lines. This value can be com- 
pared the value of 8.66 derived from the [O i] 630.03 nm line and 
the value of 8.62 derived from Oi 615.82 nm line (the latter de- 
pendent on the assumed strength of the molecular lines). The 
good performance of the 3D model in reproducing not only the 
observed line shapes and shifts, but also the centre-to-limb vari- 
ation of the oxygen lines is a good indicator of the quality of the 
inferred abundances. 



Acknowledgements. We would like to thank N. Grevesse and J. Sauval for fruit- 
ful discussions and for the use of the molecular line database. TMDP acknowl- 
edges financial support from Funda^ao para a Ciencia e Tecnologia (reference 
number SFRH/BD/2 1888/2005) and from the USO-SP International Graduate 
School for Solar Physics under a Marie Curie Early Stage Training Fellowship 
(project MEST-CT- 2005-020395) from the European Commission. This re- 
search has been partly funded by a grant from the Australian Research Council 
(DP0558836). The Swedish 1-m Solar Telescope is operated on the island of 
La Palma by the Institute for Solar Physics of the Royal Swedish Academy 



Pereira et al.: Oxygen lines in solar granulation. 11. 



15 



of Sciences in the Spanish Observatorio del Roque de los Muchachos of the 
Institute de Astroffsica de Canarias. 



References 

AUende Prieto, C, Asplund, M., & Fabiani Bendicho, P. 2004, A&A, 423, 1 109 
(AAF04) 

Allende Prieto, C, Lambert, D. L., & Asplund, M. 2001, ApJ, 556, L63 (ALAOl) 
Altrock, R. C. 1968, Sol. Phys., 5, 260 

Anders, E. & Grevesse, N. 1989, Geochim. Cosmochim. Acta, 53, 197 
Anstee, S. D. & O'Mara, B. J. 1995, MNRAS, 276, 859 
Asplund, M., Carlsson, M., & Botnen, A. V. 2003, A&A, 399, L31 
Asplund, M., Grevesse, N., Sauval, A. J., Allende Prieto, C., & Kiselman. D. 

2004, A&A, 417, 751 
Asplund, M., Grevesse, N., Sauval, J., & Scott, P 2009, ARA&A, 47, 481 
Asplund, M., Nordlund, A., Trampedach, R., Allende Prieto, C, & Stein, R. F. 

2000, A&A, 359, 729 
Ayres, T. R. 2008, ApJ, 686, 731 
Barklem, R S. 2007, A&A, 462, 781 

Barklem, R S., Belyaev, A. K., & Asplund, M. 2003, A&A, 409, LI 

Barklem, R S. & O'Mara, B. J. 1997, MNRAS, 290, 102 

Barklem, R S., O'Mara, B. J., & Ross, J. E. 1998, MNRAS, 296, 1057 

Botnen, A. 1997, Master's thesis, Inst. Theor Astrophys. Oslo (1997) 

Brault, J. W. & Neckel, H. 1987, |f tp : //f tp -hs . un i-hainburg . d e/pub/| 

outgoing/FTS-Atlas ^ 
Byrd, R. H., Lu, R, & Nocedal, J. 1995, SIAM Journal on Scientific and 

Statistical Computing, 16, 1190 
Caffau, E., Ludwig, H.-G., Stefi'en, M., et al. 2008, A&A, 488, 1031 
Centeno, R. & Socas-Navan'o, H. 2008, ApJ, 682, L61 

Cunto, W., Mendoza, C, Ochsenbein, F, & Zeippen, C. J. 1993, A&A, 275, L5+ 
Drawin, H. W. 1968, Zeitschrift fur Physik, 21 1, 404 
Eriksson, K. & Toft, S. C. 1979, A&A, 71, 178 
Eriksson, K. B. S. 1965, Ai'k. Fys., 30, 199 

Gustafsson, B., Edvardsson. B., Eriksson, K., et al. 2008, A&A, 486, 951 

Holweger, H. & Miiller, E. A. 1974, Sol. Phys., 39, 19 

Hummer, D. G., Berrington, K. A., Eissner, W., et al. 1993, A&A, 279, 298 

Johansson, S., Litzen. U., Lundberg, H., & Zhang, Z. 2003, ApJ, 584, L107 

Kiselman, D. 1993, A&A, 275, 269 

Kiselman, D. & Nordlund, A. 1995, A&A, 302, 578 

Kupka, F. G., Ryabchikova, T. A., Piskunov, N. E., Stempels, H. C., & Weiss, 

W. W. 2000, Baltic Astronomy, 9, 590 
Lambert, D. L. 1978, MNRAS, 182, 249 
Lambert, D. L. 1993, Physica Scripta Volume T, 47, 186 
Melendez, J. & Asplund, M. 2008, A&A, 490, 817 
Miiller, E. A., Baschek, B., & Holweger, H. 1968, Sol. Phys., 3, 125 
Nave, G., Johansson, S., Learner, R. C. M., Thome, A. R, & Brault. J. W. 1994, 

ApJS, 94, 221 

Nissen, R E., Primas, F, Asplund, M., & Lambert, D. L. 2002, A&A, 390, 235 
Pereira, T. M. D., Kiselman, D., & Asplund, M. 2009, A&A, in press (Paper I) 
Piskunov, N. E., Kupka, F, Ryabchikova, T. A., Weiss, W. W., & Jefl'ery, C. S. 

1995, A&AS, 112,525 
Ralchenko, Y., Kramida, A., Reader, J., & NIST ASD Team. 2009, NIST Atomic 

Spectra Database (version 3.1.5). http://physics.nist.gov/asd3 
Rosman, K. J. R. & Taylor, R D. R 1998, Pure Appl. Chem., 70, 217 
Socas-Navan-o, H. & Norton, A. A. 2007, ApJ, 660, L153 
Steenbock, W. & Holweger, H. 1984, A&A, 130, 319 
Storey, R J. & Zeippen, C. J. 2000, MNRAS, 312, 813 

Unsold, A. 1955, Physik der Stematmospharen, MIT besonderer 
Berucksichtigung der Sonne. (Berlin, Springer, 1955. 2. Aufl.) 



